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SUMMARY 

This  report  discusses  the  results  of  a  two-dimensional,  linear- 
elastic,  finite  element  analysis  of  selected  graphite/polyimide  rail 
shear  test  specimens.  The  study  includes  the  analysis  of  mechanical 
loading  and  the  effect  of  heating  the  specimen  to  a  uniform  temperature. 
The  presence  of  specimen  free  edges,  and  their  influence  on  the  accuracy 
of  the  rail  shear  test  is  discussed.  Parameters  in  this  analysis 
include  the  length-to-width  ratio  of  the  specimen  and  the  ply  layup 
for  symmetric,  balanced  laminates.  Results  presented  include  shear  and 
normal  stress  distributions  and  the  deflection  behavior  of  various 
specimens  caused  by  the  mechanical  loading  and  elevated  temperature. 


1.0  INTRODUCTION 


The  experimental  measurement  of  the  material  properties  of  laminated 
composite  materials  is  a  difficult  task.  In  particular,  both  the  deter¬ 
mination  of  the  ultimate  strength  of  a  material  in  shear  and  the  shear 
modulus  pose  problems  because  of  the  difficulties  encountered  in  attempt¬ 
ing  to  create  regions  of  pure  shear  within  a  test  specimen.  A  number  of 
experimental  methods  exist  to  determine  the  shear  properties  of  laminated 
composites.  References  1,  2  and  3  provide  discussions  and  comparisons  of 
the  variety  of  test  methods  currently  in  existence. 

Perhaps  the  simplest  and  least  expensive  of  shear  test  procedures  is 
the  so-called  rail  shear  test.  A  sketch  of  a  typical  two-rail  arrangement 
for  this  test  is  shown  in  figure  1.  The  test  method  consists  of  sand¬ 
wiching  a  small,  flat,  rectangular  test  specimen  between  large  metallic 
rails.  The  attachment  of  the  test  specimen  to  the  individual  rails  is 
accomplished  either  by  a  clamping  action  provided  by  bolting  the  rails 
to  the  specimen  or  by  adhesive  bonding  between  the  rails  and  the  specimen. 
The  load  from  the  test  machine  can  be  introduced  into  the  specimen  in 
one  of  two  ways.  One  method  of  load  introduction  is  to  apply  tensile 
or  compressive  loads  along  the  specimen  centerline  parallel  to  the  rails 
themselves,  as  indicated  in  figure  1.  An  inplane  shear  load  is  trans¬ 
mitted  to  the  test  specimen  by  displacing  one  rail  relative  to  the  other. 
An  alternative  method  of  load  application  is  to  apply  the  loads  along  a 
diagonal  of  the  rectangular  specimen. 

In  any  test  to  determine  either  shear  modulus  or  shear  strength  of 
a  laminated  composite  material,  a  state  of  uniform  shear  stress  should 
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exist  within  the  boundaries  of  the  test  specimen.  For  a  flat  rectangu¬ 
lar  specimen,  a  uniform  shear  state  can  be  imposed  only  by  loading  all 
four  boundaries;  such  simultaneous  loading  of  the  specimen  boundaries  is 
not  present  in  the  standard  rail  shear  test.  A  result  of  the  rail  shear 
type  of  loading  is  that  two  parallel  edges  of  the  specimen  are  loaded 
while  the  remaining  two  parallel  edges  remain  stress  free.  It  has  been 
shown  (ref.  4)  that  the  presence  of  these  free  edges  gives  rise  to  edge 
zones  in  which  the  inplane  shear  stress  varies  greatly.  Accompanying 
these  shear  stress  variations  are  normal  stresses  that  ensure  stress 
equilibrium.  The  fact  that  normal  stresses  are  developed  within  the 
specimen  by  the  loading  shown  in  figure  1  has  given  rise  to  criticism 
of  the  use  of  the  rail  shear  test  to  determine  laminate  properties 
(ref.  5).  These  normal  stresses  may  cause  undesirable  occurrences  such 
as  debonding  of  the  specimen/rail  interface  or  failure  of  the  specimen 
near  its  corners  due  to  mechanisms  other  than  shear  initiated  failure. 

In  addition,  the  presence  of  a  nonuniform  shear  stress  distribution 
within  the  specimen  may  lead  to  considerable  error  in  the  determination 
of  the  shear  modulus. 

Previous  studies,  such  as  those  by  Whitney,  Stansbarger  and 
Howell  (ref.  5)  and  Whitney  (refs.  6,7),  have  shown  that  the  nonuniformity 
of  the  stresses  in  the  edge  zone,  and  thus  the  deviation  from  a  state  of 
pure  shear  within  the  specimen,  is  governed  by  such  factors  as  the  speci¬ 
men  length-to-width  ratio  (also  referred  to  as  the  aspect  ratio,  b/a  in 
figure  2)  and  the  laminate  layup. 

An  additional  problem  appears  when  testing  specimens  at  elevated 
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temperatures.  Because  the  rails  and  the  composite  specimen  are  made  of 
dissimilar  materials,  thermally  induced  stresses  are  caused  by  differential 
thermal  expansion.  This  problem  is  discussed  by  Bergner,  Davis  and 
Herakovich  (ref.  3).  Results  of  a  finite  element  analysis  of  a  large 
length-to-width  ratio  specimen  are  presented  in  this  latter  report; 
these  results  show  that  elevated  temperatures  may  produce  additional 
large  normal  stresses  near  the  test  specimen  corners.  The  analytical 
results  presented  in  reference  3  are  limited  in  scope  since  the  effect 
of  specimen  aspect  ratio  on  test  validity  is  not  considered. 

In  an  effort  to  obtain  more  comprehensive  knowledge  and  under¬ 
standing  of  the  stress  fields  that  arise  in  a  rail  shear  test  with 
loading  of  the  type  illustrated  in  figure  1,  both  at  room  temperature 
and  at  elevated  temperature,  a  series  of  analytical  studies  using  finite 
element  models  was  performed.  The  analyses  and  results  presented  in 
this  report  are  those  that  consider  a  two-dimensional,  balanced, 
symmetrical -ply,  rail  shear  test  specimen  with  flexible  rails.  Classical 
linear-elastic,  laminated  plate  theory  is  used  for  the  finite  element 
model.  Specimen  aspect  ratios  of  between  2  and  12  are  considered. 

Results  presented  include  shear  and  normal  stress  distribution  infor¬ 
mation  at  various  positions  on  the  rail  shear  specimen,  both  for  mech¬ 
anical  and  thermal  loading.  For  the  mechanical  loading  case,  the 
laminates  considered  are  [0^2^’  [90i2]>  [90/145/0]^  and  [±462]^ 
layups,  while  the  thermal  load  cases  include  only  the  [O12]  and  [9O12] 
laminates.  From  these  results,  conclusions  are  drawn  about  the 
acceptability  of  the  rail  shear  test  as  an  experimental  tool  for  shear 
testing  of  composites. 
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2.0  SYMBOLS 


Definitions  of  the  major  symbols  used  in  this  report  are  given  in 
this  section.  Symbols  not  having  general  use  are  defined  as  they  are 
introduced.  The  units  used  for  physical  quantities  are  given  in  the 
International  System  of  Units  (SI),  except  where  noted. 

a,b  Specimen  width  and  length  respectively 

(figure  2). 

Eo  Modulus  of  elasticity  in  the  laminate  0° 

direction. 

P  Test  machine  load  applied  to  rail  (figure  1). 

t  Test  specimen  thickness. 

AT  Change  in  temperature. 

Coefficient  of  thermal  expansion  of  the  rail. 

a  Coefficient  of  thermal  expansion  of  the 

^  specimen  in  the  0°  direction. 

T  Shear  stress  in  the  x-y  coordinate  system. 

T  Average  shear  stress,  x  =  P/bt. 

Oy  Reference  stress  for  thermal  stress 

results,  Eqn.  1. 
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3.0  IDEALIZATION  AND  ANALYSIS  OF 
RAIL  SHEAR  TEST  SPECIMENS 

3. 1  Background  -  Analysis  of  Mechanical  Loading 

The  interest  in  this  report  is  focused  upon  the  two-dimensional 
stress  distributions  in  flat  rectangular  specimens  that  have  balanced, 
symmetrical -ply,  laminated  composite  material  and  are  subjected  to  in- 
plane  shear  loading.  Furthermore,  the  mechanical  loading  may  occur 
both  at  room  temperature  and  at  elevated  temperature.  Because  the 
analytical  model  to  be  discussed  is  linear  elastic,  superposition  prin¬ 
ciples  apply.  Thus  the  overall  problem  may  be  subdivided  into  two 
parts:  (1)  mechanical  loading  with  no  temperature  effects  included; 
and  (2)  the  introduction  of  temperature  effects  with  no  external  load. 

The  interest  in  the  mechanical  loading  of  flat  rectangular  plates  by 
inplane  shear  loads  of  the  type  found  in  rail  shear  testing  extends  back 
to  the  early  part  of  this  century.  In  1912,  Coker  (ref.  4)  conducted 
photoelastic  experiments  to  determine  the  effect  of  specimen  aspect  ratio 
(length-to-width  ratio)  on  the  shear  stress  distribution  in  isotropic 
flat  plates  with  parallel  edges  subjected  to  edge  shear,  a  loading  that 
is  nearly  identical  to  that  imposed  during  the  rail  shear  test.  For  small 
aspect  ratio  specimens,  Coker  determined  a  shear  stress  distribution, 
along  the  specimen  centerline  parallel  to  the  applied  shear  load,  much 
like  the  parabolic  distribution  that  is  predicted  from  elementary, 
strength  of  material,  beam  theory.  This  parabolic  distribution  of  shear 
stress  was  found  to  change  as  the  specimen  aspect  ratio  increased.  In¬ 
stead  of  a  single  maximum  stress,  as  occurs  for  a  parabolic  shear  stress 
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distribution,  plates  with  larger  aspect  ratios  (of  the  order  of  4) 
exhibit  two  shear  stress  maxima,  symmetrically  located  about  the  plate 
center.  These  maxima  are  located  at  a  distance  of  about  3/4  of  an 
edge  width  from  the  free  edges  of  the  plate  specimen.  Coker  further 
showed,  by  photoelastic  experiment,  that  the  shear  stress  remains  at 
a  reasonably  constant  value  in  the  interior  of  the  plate.  This  latter 
result  of  Coker's  work  is  frequently  cited  as  justification  for  the  use 
of  large  aspect  ratio,  rail  shear  specimens  for  testing  of  isotropic 
materials. 

In  1923,  Inglis  (ref.  8)  published  an  analytical  study  of  Coker's 
experimental  work.  In  this  study,  Inglis  used  a  classical  stress  function 
approach  to  formulate  and  to  solve,  approximately,  the  theory  of 
elasticity  problem  for  a  linear  elastic,  isotropic  material  with  two 
parallel  edges  loaded  in  shear  and  the  other  two  edges  free.  The  loaded 
boundaries  are  assumed  to  be  rigid.  A  series  solution  was  obtained  that 
is  quite  accurate  in  regions  away  from  the  plate  corners.  Inglis  found 
that  a  large  stress,  normal  to  the  boundary  loaded  in  shear,  appears  near 
the  corners  of  the  plate.  This  normal  stress  has  a  magnitude  four  times 
that  of  the  average  inplane  shear  stress,  x,  applied  to  the  plate 
boundaries.  It  was  noted  that  this  stress  concentration  factor  of  four 
is  somewhat  inaccurate  because  of  the  truncation  of  the  assumed  series 
solution. 

The  analysis  of  two  other  problems  related  to  the  rail  shear  test 
received  considerable  attention  in  the  late  1930's  and  during  the  1940's. 
The  first  problem,  that  of  shear  lag  in  reinforced  sheet-stringer 


9 


construction,  was  the  subject  of  a  number  of  papers  and  reports  (see  Kuhn, 
ref.  9).  The  shear  lag  problem  approximates  the  theory  of  elasticity 
problem  in  which  axial  loads  are  transferred  from  flexible  longitudinal 
stiffeners  to  thin  sheets  of  shear  web  material  attached  to  the  stiffeners. 
The  key  assumption  of  shear  lag  analysis  is  that  the  transverse  stiffness 
of  the  specimen  is  infinite.  The  important  finding  of  these  shear  lag 
analyses  for  isotropic  materials  is  that  the  load  transfer,  or  load  dif¬ 
fusion,  from  flexible  stiffeners  to  the  shear  web  material  produces 
regions  of  non-uniform  stress  at  the  panel  ends  where  the  axial  load  is 
introduced.  This  non-uniform  stress  distribution  decays  exponentially 
with  distance  if  the  panel  is  long  enough;  the  argument  of  the  exponential 
decay  function  is  related  to  stiffener  and  sheet  material  properties  and 

1  /2 

sheet  geometry.  In  particular,  for  isotropic  materials  the  factor  (G/E) 
is  found  to  be  a  parameter  in  this  exponential  argument.  The  parameter  G 
is  the  shear  modulus  of  the  material  while  E  represents  the  modulus  of 
elasticity  in  the  direction  of  the  applied  axial  load. 

Hildebrand  (ref.  10)  presents  several  exact  solutions  to  the  differ¬ 
ential  equations  that  govern  shear-lag  in  orthotropic  panels  that  are 
assumed  to  be  rigid  in  the  transverse  direction  (perpendicular  to  the 
applied  axial  load  direction).  While  this  latter  assumption  precludes 
the  use  of  these  results  for  the  rail  shear  problem,  Hildebrand's  results 
show  the  exponential  decay  character  of  the  stresses  caused  by  this  type 
of  loading. 

A  field  of  study  with  characteristics  related  to  the  rail  shear 
test  is  found  in  the  area  of  bonded  joints.  In  particular,  the  work  of 
Goland  and  Reissner  (ref.  11)  predicts  normal  stresses  at  corners  of 
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bonded  lap  joints.  These  "tearing"  stresses  are  set  up  in  the  adhesive 
as  a  result  of  the  eccentricity  of  the  loading  on  the  lap  joint.  Goland 
and  Reissner  make  certain  simplifying  assumptions  which  cause  their 
results  to  be  not  strictly  applicable  for  analysis  of  the  rail  shear 
test.  For  instance,  the  adherends  are  assumed  to  be  relatively  flexible 
and  adhesive  flexibility  is  ignored  in  a  direction  transverse  to  the 
load  direction.  Their  studies  do,  however,  reveal  that  large  tearing 
stresses  normal  to  the  applied  load  direction  may  be  expected  in  the 
vicinity  of  the  joint  edges.  In  a  rail  shear  test,  these  joint  edges 
correspond  to  the  rail  shear  specimen  corners. 

More  recently,  the  application  of  the  rail -shear  test  to  advanced 
composite  materials  has  received  considerable  attention.  Whitney, 
Stansbarger  and  Howell  (ref.  5)  and  Whitney  (references  6,  7)  have 
studied  the  rail  shear  test  in  some  detail.  Reference  6  presents  a 
detailed  analysis  of  the  stresses  in  a  symmetrical,  balanced-ply, 
laminated  plate  loaded  by  inplane  shear  along  two  parallel  edges  that 
are  rigid  but  displace  relative  to  each  other  to  cause  shear  in  the 
specimen.  One  potential  difficulty  with  the  analysis  presented  is  that  a 
Fourier  series  solution  for  the  displacement  field  was  used  to  solve  the 
problem.  The  series  solution  used  is  discontinuous  at  the  free  edges; 
thus,  the  solution  discontinuity  at  the  edge  may  be  masking  the  free 
edge  effects.  In  reference  7,  test  results  are  presented  for  a  rec¬ 
tangular  orthotropic  rail  shear  specimen  clamped  to  metallic  rails, 
while  reference  5  summarizes  much  of  the  work  presented  in  references  6 
and  7. 

Bergner,  Davis  and  Herakovich  (ref.  3)  present  a  limited  analytical 
study  of  the  rail  shear  test  with  composite  materials.  Reference  3  is 
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significant  because,  unlike  the  differential  equation  approach  used  in 
previous  studies,  the  authors  use  a  finite  element  analysis  for  the 
study  of  rail  shear  testing  both  at  room  temperature  and  also  with 
temperature  effects.  Their  results  show  the  superiority  of  the  finite 
element  approach  over  the  Fourier  series  solution. 

3.2  Background  -  Analysis  of  Elevated  Temperature  Tests 

The  elevated  temperature  rail  shear  test  has  not  been  examined  in 
great  detail.  Reference  3  presents  an  analysis  of  the  stresses  occurring 
in  a  rail  shear  test  specimen  with  flexible  rails  when  the  temperature 
is  lowered  316°C  (600°F).  The  results  of  that  study  show  that  regions 
of  large  normal  stress  occur  in  the  specimen  corners  as  the  tempera¬ 
ture  is  changed;  this  thermal  stress  is  found  to  be  significant. 

Two  previous  studies  consider  a  problem  that  is  similar  to  that 
encountered  in  elevated  temperature  rail  shear  testing.  Aleck  (ref.  12) 
has  analyzed  the  thermal  stresses  in  a  rectangular  plate  of  isotropic 
material  in  which  one  edge  is  cemented  to  a  perfectly  rigid  plane  while 
the  other  three  edges  are  allowed  to  be  free.  The  solution  of  this 
problem  with  energy  methods  yielded  a  solution  that  is  suspect,  in  that 
the  stress  functions  used  to  approximate  the  solution  are  somewhat 
simplified  and  thus  inadequate  to  represent  the  solution. 

A  thermal  stress  study  by  Kobatake  and  Inoue  (ref.  13)  concerns 
a  linear  elastic,  isotropic  rectangular  plate,  two  of  whose  parallel 
edges  have  been  fixed  to  rigid  planes  before  heating  is  introduced. 

The  study  presented  in  reference  13  indicates  that,  due  to  the  restrained 
thermal  expansion,  large  normal  stresses  and  shear  stresses  appear  near 
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the  plate  corners.  However,  the  tnagnitude  of  these  stresses  tends  to 
decrease  rapidly  with  distance  away  fron  the  free  edges  (distances  that 
are  slightly  greater  than  an  edge  width  of  the  plate).  Unfortunately, 
boundary  conditions  are  applied  to  the  plate  that  lead  to  edges  free 
of  stress  normal  to  the  interface.  As  a  result  of  this,  the  stress 
distributions  near  the  interface,  presented  in  reference  13,  are  not 
correct  for  the  present  problem.  Nevertheless,  reference  13  indicates 
that  sizable  shear  and  normal  stresses  near  the  specimen  corners  may 
be  caused  by  constrained  expansion  of  the  rails  and  the  specimen  during 
an  elevated  temperature  rail  shear  test.  For  instance,  the  maximum 
value  of  shear  stress  computed  near  the  specimen  corners  is,  for  their 
isotropic  material,  x  =  0.645  E(Aa)  (aT),  where  E,  Aa  and  AT  represent 
the  modulus  of  elasticity,  the  difference  between  the  values  of  the 
coefficients  of  thermal  expansion  of  the  rails  and  specimen,  and  the 
change  in  temperature,  respectively.  Of  further  significance,  is  that 
the  shear  stress  caused  by  thermal  expansion  is  minimal  at  a  distance 
of  two  plate  edge  widths  into  the  plate  from  the  free  edges  and,  at 
one  width,  has  a  value  of  only  about  25%  of  the  maximum  value  given 
above.  Symmetry  considerations  indicate  that  the  shear  stress  x  is 
zero  along  both  specimen  axes  of  synmetry. 

A  conclusion  to  be  drawn  from  reference  13  is  that,  during  a 
rail  shear  test  at  elevated  temperature,  stresses  caused  by  thermal 
expansion  probably  will  not  affect  significantly  the  determination  of 
the  shear  modulus.  However,  the  shear  and  normal  stresses  generated 
near  the  corners  of  the  specimen  at  the  specimen/rail  interface  may 
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lead  to  premature  failure  of  the  specimen  at  or  near  this  interface. 

For  adhesively  bonded,  rail /shear  specimens,  these  stresses  may 
cause  the  bond  itself  to  fail  due  to  insufficient  strength. 

The  studies  previously  cited  provide  guidance  for  the  idealization 
and  stress  analysis  of  the  rail  shear  test  specimen.  These  studies 
show  that  a  finite  element  model  provides  a  simpler  and  more  accurate 
analytical  tool  for  the  study.  In  addition,  the  flexibility  of  the 
rails  should  be  considered.  Finally,  careful  attention  must  be  given  to 
the  model  idealization  in  the  vicinity  of  the  specimen  corners  where 
large  stress  gradients  will  occur.  The  next  section  describes  the 
finite  element  model  used  for  this  study. 

3.3  Model  Idealization  and  Analysis 

The  formulation  of  an  accurate  but  economical  (in  terms  of  computer 
time)  finite  element  model  of  a  rail/specimen  combination,  such  as  shown 
in  figure  1,  requires  that  some  assumptions  be  made  about  the  expected 
physical  behavior  of  the  test  specimen  and  its  rails  under  the  applied 
load.  The  idealized  model  chosen  to  approximate  this  behavior  focuses 
attention  on  the  two-dimensional  displacement  and  stress  behavior  of  the 
composite  test  specimen,  while  retaining  the  important  influence  of  the 
axial  flexibility  and  bending  flexibility  of  the  two  titanium  rails. 

Consider  figure  3,  in  which  is  shown  a  typical  finite  element 
idealization  of  the  rail  shear  test  configuration  used  in  this  analysis. 
The  flat  rectangular  specimen  itself  is  modelled  as  an  assemblage  of 
planar,  rectangular  finite  elements  whose  material  behavior  is  linear 
elastic,  but  anisotropic.  The  rails  are  titanium,  chosen  because  of 
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Figure  3  -  Finite  Element  Model  of  Rail/Specimen  Combination. 

(Note  that  x-y  coordinate  origin  is  in  the  lower 
left  corner  of  the  specimen.) 


their  characteristic  small  coefficient  of  thermal  expansion.  These 
rails  are  represented  by  beam  elements  with  a  straight  elastic  axis  in 
the  plane  of  the  specimen  and  parallel  to  the  y-axis  of  the  specimen  in 
figure  3. 

Although  the  load  from  the  test  machine  is  applied  along  the 
specimen  centerline  parallel  to  the  y-axis  in  figure  3,  the  mathematical 
idealization  employs  a  statically  equivalent  load  combination  consisting 
of  a  concentrated  force  and  moment  applied  to  the  end  of  a  straight 
beam  that  simulates  the  rail.  This  load  set  is  illustrated  in  figure  2. 

Since  the  finite  element  analysis  involves  a  numerical  representa¬ 
tion  of  the  structure,  while  the  interpretation  of  the  analysis  results 
is  best  accomplished  in  nondimensional  parameter  form,  consistent 
parameters  are  assigned  to  the  models  so  that  some  parameters  are 
fixed  while  others  are  varied.  From  the  literature  previously  cited,  it 
is  expected  that  two  effects  will  affect  the  diffusion  of  the  axial 
loads  from  the  rails  into  the  test  specimen.  These  effects  are  the 
presence  of  the  specimen  free  edges  and  the  elastic  interaction  between 
the  rails  and  the  specimen.  For  these  reasons,  the  cross-sectional  area 
and  shape  of  the  rails  was  held  fixed,  as  was  the  specimen  width,  the 
dimension  "a"  in  figure  2.  For  aspect  ratio  studies,  the  dimension  b 
was  varied  to  achieve  the  proper  b/a  ratio. 

Idealized  structural  models  with  geometries  of  the  type  shown  in 
figure  3  were  developed  and  analyzed  using  the  SPAR  computer  program 
(ref.  14)  at  the  NASA/Langley  Research  Center.  The  results  of  these 
analyses  are  described  in  Section  4.0.  The  geometrical  parameters  and 
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material  properties  of  the  rails  and  the  graphite/polyimide  material 
are  given  in  the  Appendix. 
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4.0  RESULTS  AND  DISCUSSION 


The  importance  of  the  rail  shear  test  specimen  aspect  ratio  and  the 
effect  of  the  laminate  layup  on  the  stress  field  present  during  testing 
have  been  noted  in  the  previous  discussion.  To  study  the  effects  of  aspect 
ratio  and  laminate  layup,  several  different  aspect  ratios  were  considered 
for  four  types  of  laminates.  The  laminates  are  composed  of  multiple 
layers  of  graphite/polyimide  material.  The  laminates  considered  are: 

(1)  a  ^^-[2} 

(2)  a  [90^ 2^  laminate, 

(3)  an  8-layer  [±45215  laminate, 

(4)  an  8-layer  [90/±45/0]^  laminate. 

The  reference  position  for  the  measurement  of  the  fiber  angle  e  is 
illustrated  in  figure  2. 

4. 1  Mechanical  Loading 

Considering  first  a  specimen  with  an  aspect  ratio  of  2,  it  is  seen 
in  figure  4  that  laminate  construction  greatly  affects  the  shear  stress 
distribution  along  the  centerline  parallel  to  the  y-axis.  Figure  4 
displays  the  laminate  shear  stress  versus  distance  from  the  free  edge; 
this  latter  distance  is  nondimensional ized  with  respect  to  the  specimen 
width.  The  shear  stress  in  the  laminate  is  nondimensional ized  with 
respect  to  the  average  value  of  the  shear  stress,  computed  as  t  =  P/bt. 
Since  these  distributions  of  shear  stress  are  symmetrical  about  y  =  b/2, 
only  the  distributions  present  in  half  the  specimen  are  shown. 

It  must  be  noted  here  that  the  stresses  shown  throughout  this 
report  are  the  result  of  taking  the  resultant  inplane  forces  per 
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Figure  4  -  Shear  Stress  Ratio  t/t  versus  y/a  at  x=a/2  for 
Four  Different  Specimens.  Aspect  Ratio  Equals 


unit  length,  N  ,  N  and  present  in  the  laminate  and  dividing  by  the 
laminate  thickness  t.  Thus,  depending  upon  the  laminate  layup,  the 
stresses  in  the  individual  lamina  may  be  different  than  those  values 
given  in  the  ensuing  discussion. 

For  the  [90.|2]  laminate,  the  discrepancy  between  the  actual  value 
of  shear  stress  t  and  the  average  value  of  shear  stress,  x,  is  seen  to 
be  5%  at  the  specimen  center.  The  discrepancy  is  50%  for  the  [±452]g 
laminate.  In  addition,  the  shapes  of  the  nondimensional  shear  stress 
distributions  are  different.  The  shape  of  the  x/x  distribution  for  the 
[0-J2]  laminate  and  that  for  the  [90/±45/0]^  laminate  more  nearly 
resemble  the  parabolic  distributions  observed  by  Coker  for  isotropic 
specimens  with  the  same  aspect  ratio. 

When  the  aspect  ratio  of  the  specimen  is  increased  to  3,  the  disparity 
between  the  true  shear  stress  and  the  average  shear  stress  at  the  specimen 
center  is  less  than  that  previously  seen,  as  is  shown  in  figure  5.  The 
discrepancy  at  the  specimen  center  between  the  average  shear  stress  and 
the  true  shear  stress  is  as  little  as  3%  for  the  [9O12]  laminate.  The 
r±45o]  laminate  still  shows  a  27%  difference  at  the  center.  The  data 
for  the  [90/±45/0]^  laminate  is  not  shown  because  its  behavior  so 
closely  approximates  that  of  the  [O12]  laminate. 

Figure  6  illustrates  the  effect  on  the  shear  stress  distribution 
of  increasing  still  further  the  specimen  aspect  ratio  to  4.  Differences 
between  the  average  shear  stress  x  and  the  true  shear  stress  x  at  the 
specimen  center  decrease  from  those  shown  in  figure  5  and  are  between 
2%  and  17%.  In  addition,  this  figure  clearly  shows  the  presence  of  a 
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.27,  [±45,]  Symmetrical 


Laminate 


local  rnaxiitium  shear  stress,  a  type  of  shear  stress  distribution  similar 
to  that  noted  by  Coker  in  his  results  for  specimens  with  a  similar 
aspect  ratio. 

By  further  increasing  the  aspect  ratio  from  4  to  a  value  of  6,  the 
difference  between  the  average  shear  stress  and  the  true  shear  stress 
narrows  even  more.  Figure  7  shows  these  differences  at  the  specimen 
center  to  be  from  between  ^%  and  8%,  depending  upon  the  type  of  laminate. 

To  illustrate  further  the  effect  of  changing  specimen  aspect  ratio 
upon  the  shear  stress  distribution  within  the  specimen,  let  us  focus 
attention  upon  two  laminates,  the  [90^2^  laminate  and  the  [±452]^ 
laminate.  These  laminates  are  chosen  because  of  the  contrasting  be¬ 
havior  of  their  shear  stress  distributions. 

Figures  8  and  9  illustrate  the  behavior  of  the  shear  stress 
distribution  for  these  two  laminates  as  specimen  aspect  ratio  is 
increased  from  6  to  12.  In  figure  8,  it  is  seen  that  doubling  the 
aspect  ratio  of  the  [9O12]  specimen  from  6  to  12  changes  very  little 
the  shear  stress  distribution  along  the  specimen  centerline.  The  shear 
stress  ratio  t/t  at  the  specimen  center  is  less  than  unity  for  an 
aspect  ratio  of  12.  As  a  consequence,  the  maximum  value  of  the  shear 
stress,  found  near  the  free  edge,  increases  as  the  aspect  ratio  in¬ 
creases.  However,  the  change  in  these  maximum  values  is  slight. 

Figure  9  shows  shear  stress  distributions  for  the  [±^52]^  laminate 
that  are  in  marked  contrast  to  those  seen  in  figure  8.  As  specimen 
aspect  ratio  increases,  a  very  pronounced  increase  in  the  maximum  shear 
stress,  occurring  near  the  free  edge,  is  seen.  Even  so,  the  shear 
stress  ratio  t/t,  present  at  the  specimen  center,  is  very  nearly  unity. 
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Figure  7  -  Shear  Stress  Ratio  t/t  versus  y/a  at  x=a/2  for 
Three  Different  Specimens.  Aspect  Ratio  Equals 


versus  y/a  at  x=a/2;  [90 


ure  9  -  Shear  Stress  Ratio  t/t  versus  y/a  at  x=a 
Laminate,  Aspect  Ratios  6,8,12. 


The  appearance  of  a  nonuniform  shear  stress  field  within  the 
specimen  gives  rise  to  normal  stress  fields.  An  example  of  the  a  stress 
distribution  that  occurs  is  shown  in  figure  10.  In  figure  10,  data 
taken  from  the  finite  element  analysis  of  a  [90^ 2^  laminate  with  aspect 
ratio  6  is  displayed.  This  curve  represents  the  normal  stress  a 

A 

(normalized  with  respect  to  the  average  shear  stress  x),  occurring  at 
the  centers  of  finite  elements  used  in  the  specimen  model,  adjacent  to 
the  rail/specimen  interface.  For  this  model,  the  locus  of  these  points 
is  a  line  at  a  distance  0.025a  from  this  interface.  The  stress 
distribution  shown  in  figure  10  illustrates  the  large  stress  concentra¬ 
tions  present  at  the  specimen  corners.  The  normal  stress  is  seen  to 
decrease  rapidly  with  distance  from  the  free  edge. 

The  ratio  a  /x  calculated  near  the  corners  of  specimens  is  shown 

A 

in  the  table  below  for  several  aspect  ratios. 


Location 
of  Data 
(x,y) 

1 — I 

O 

ro 

1 _ t 

CO 

1 - \ 

OJ 

LO 

+1 

1 - 1 

[O12] 

[90/±45/0]^ 

2 

(0.025a, 0.025a) 

6.97 

2.31 

3.05 

3.35 

3 

(0.025a, 0.025a) 

6.95 

2.03 

2.88 

3.21 

4 

(0.025a, 0.067a) 

6.14 

2.63 

2.97 

6 

(0.025a, 0.025a) 

7.51 

2.07 

2.83 

3.40 

10 

(0.05a, 0.025a) 

8.46 

2.44 

2.73 

3.97 

STRESS  CONCENTRATION  FACTORS  o/x  NEAR  SPECIMEN  CORNERS 

A 


27 


29 


Figure  11  -  Mechanically  Loaded  Specimen  Deformation  Figure  12  -  Mechanically  Loaded  Specimen  Deformation 
Pattern;  [DO-jg]  Laminate;  Aspect  Ratio  Pattern;  [±4523s  Laminate;  Aspect  Ratio 

Equals  2.  Each  Grid  Section  Encloses  Equals  2.  Each  Grid  Section  Encloses 

Four  Finite  Elements.  Four  Finite  Elements. 
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Deformation  Pattern;  [9O12]  Deformation  Pattern;  [±45215 

Laminate;  Aspect  Ratio  Equals  Laminate;  Aspect  Ratio  Equals 
4.  Each  Grid  Section  Encloses  4.  Each  Grid  Section  Encloses 
Four  Finite  Elements.  Four  Finite  Elements. 


From  this  table  of  values  it  is  seen  that  the  ratio  o  /x  near  a 

X 

specimen  corner  can  be  quite  large,  particularly  in  the  case  of  the 
[9O12]  specimen.  Note  that  the  data  displayed  in  the  table  are  not 
given  at  exactly  the  same  x,y  positions  due  to  differences  in  finite 
element  mesh  sizes  for  the  various  aspect  ratios. 

To  complete  the  presentation  of  the  effects  of  mechanical  loading 
upon  the  rail  shear  specimen,  figures  11  through  14  are  presented; 
these  figures  contrast  the  displacement  behavior  of  the  [90.J2]  laminate 
with  the  [±462]^  laminate  for  aspect  ratios  of  2  and  4. 

4.2  Differential  Thermal  Expansion 

Because  of  the  mismatch  between  coefficients  of  thermal  expansion, 
the  heating  or  cooling  of  a  rail  shear  specimen  will  cause  stresses  to 
develop  in  the  rails  and  the  test  specimen.  Fortunately,  the  geometrical 
symmetry  of  the  specimen  and  the  rails  and  the  nature  of  the  thermal 
loading  preclude  the  development  of  a  shear  stress  x  along  either  of  the 
centerlines  of  the  specimen.  For  this  reason,  the  experimental  deter¬ 
mination  of  the  shear  modulus  is  unlikely  to  be  affected  by  temperature 
if  the  shear  strain  is  measured  at  the  specimen  center.  However,  near 
the  corners  of  the  specimen,  sizable  thermal  stresses  may  develop, 
stresses  that  may  cause  experimental  error  when  testing  for  the  ultimate 
strength  of  the  laminate  in  shear. 

The  finite  element  models  used  in  this  portion  of  the  study  assume 
the  rail  and  specimen  to  be  heated  to  a  uniform  temperature.  In 
addition,  the  specimen  and  the  attached  rails  are  free  to  expand 
laterally.  The  finite  element  idealizations  used  for  the  thermal 
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expansion  study  are  identical  to  those  described  in  the  previous  study 
of  mechanical  loading. 

To  illustrate  the  response  of  an  orthotropic  test  specimen  to  a 
uniform  temperature  increase,  results  are  presented  for  two  test  speci¬ 
mens,  the  [0^2^  [90^2^  laminates;  four  aspect  ratios  are  analyzed. 

Nondimensional ization  of  the  results  is  accomplished  by  dividing  the 
thermally  induced  stress  at  any  point  on  the  laminate  by  a  reference 
stress  Oj,  given  by  the  expression 

Oj  =  Ep  (Aa)  (aT)  ) 

In  Eqn.  1,  E^  is  the  modulus  of  elasticity  of  the  specimen  in  the  0° 
direction  of  the  laminate,  while  Aa  is  given  by  the  relationship 

Aa  =  (2) 

In  Eqn.  2,  a^  is  the  coefficient  of  thermal  expansion  of  the  rails, 

while  a  is  the  coefficient  of  thermal  expansion  of  the  specimen  in  the 
s 

0°  direction.  AT  is  defined  as  the  temperature  change  from  the  zero 
strain  condition.  Note  that  the  value  of  depends  upon  the  specimen 
being  studied. 

A  physical  interpretation  of  is  that  it  is  the  value  of  the 
compressive  stress  necessary  to  keep  the  extensional  strain  in  the  0° 
direction  of  the  specimen  equal  to  zero  for  a  material  with  a 
coefficient  of  thermal  expansion  of  Aa  and  a  temperature  change  AT. 
Since  no  such  stress  can  be  applied  to  the  specimen,  because  of  the 
two  stress  free  edges,  the  specimen  normal  stresses  near  these  edges 
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in  the  y  direction  will  differ  from  this  value,  but  it  will  be  shown  that 
the  stresses  will  approach  this  value  over  the  central  region  of 
the  specimen  if  the  aspect  ratio  is  large  enough. 

4.2.1  Results  for  the  [9O12]  Laminate 

Shown  in  figure  15  are  the  normalized  distribution?  of  stress 

and  Oy/oj  that  occur  along  the  line  x=a/2  for  a  [90.|2]  laminate  with  an 

aspect  ratio  of  2.  If  the  specimen  is  heated,  the  stress  will  attain 

a  maximum  compressive  value  near  the  center  of  the  specimen.  At  the 

center  of  the  free  edge  of  the  specimen,  the  normal  stress  a  will  have 

a  large  tensile  value,  a  value  that  declines  rapidly  with  distance  away 

from  the  free  edge.  It  should  be  noted  that  the  reference  stress  a-j- 

for  this  laminate  will  be  negative  for  an  increase  in  temperature,  due 

to  the  fact  that  exceeds  a„. 

s  r 

Figures  16  and  17  show  the  behavior  of  the  normal  stresses 
and  a^,  at  x=a/2,  for  the  [90.|2]  laminate  with  specimen  aspect  ratios 
of  4  and  10,  respectively.  These  stresses,  nondimensional ized  with 
respect  to  CTj  for  this  laminate  and  plotted  against  y/a,  behave  in 
a  manner  similar  to  those  plotted  in  figure  15. 

The  behavior  of  the  normal  stress  a  near  the  free  edge  of  a 

A 

[90^12]  specimen  with  an  aspect  ratio  of  10  is  illustrated  in  figure  18. 

In  this  figure,  the  quantity  at  the  position  y/a"0.025  is  graphed 

versus  x/a.  This  figure  shows  that,  if  the  specimen  is  heated,  a  large 
compressive  stress,  perpendicular  to  the  rails,  will  exist  near  the 
specimen  corners,  positions  corresponding  to  the  locations  x=0  and  x=a 
in  figure  18. 
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Figure  15  -  Nondimensional  Normal  Stresses  and  Oy/oj  versus 

y/a  at  x=a/2.  [9O12]  Laminate;  Specimen  Aspect  Ratio 
Equals  2. 
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Figure  16  -  Nondimensional  Normal  Stresses  a  /oy  and  o  /a^  versus  y/a  at  x=a/2. 


armnate;  Specimen 


To  illustrate  the  nonuniformity  of  the  normal  stresses  that  may 
exist  at  the  specimen  centerline  position  y=b/2,  figure  19  presents  the 
nondimensional ized  stresses  and  ^s  a  function  of  x/a  at 

y=b/2,  for  specimen  aspect  ratios  of  2,  4  and  6.  For  specimen  aspect 
ratios  greater  than  6,  the  normal  stresses  are  constant  across  the 
specimen  at  y=b/2  and  cannot  be  distinguished  from  those  given  for  an 
aspect  ratio  of  6. 

4.2.2  Results  for  the  [O-j 2]  Laminate 

If  the  specimen  fibers  are  oriented  parallel  to  the  rails,  the 
thermal  stress  problem  changes.  To  illustrate  the  effect  of  a  uniform 
at  on  a  specimen  whose  fibers  are  oriented  parallel  to  the  rails,  a 
[O12]  laminate  was  studied.  The  coefficient  of  thermal  expansion  of 
the  composite  specimen  in  the  direction  of  the  rails  is  nearly  zero 
so  that  the  rail/specimen  interface  tends  to  elongate  with  increases 
in  temperature.  This  elongation  causes  tensile  stresses  in  the 
y-direction  near  the  center  of  the  specimen.  Note  that,  for  this 
laminate,  Oj  will  be  a  positive  number  if  the  specimen  is  heated. 

Figure  20  shows  the  normal  stresses  that  occur  along  a  line  at 
x=a/2  for  a  [0^2^  specimen  with  aspect  ratio  of  2.  When  compared  to 
figure  15,  this  figure  shows  that  the  stresses  in  the  [0.|2]  specimen 
decrease  less  rapidly  with  distance  away  from  the  free  edge.  In  addition, 
the  cr^  stress  distribution  shown  in  figure  20  changes  sign,  unlike  the 
CTy  distribution  shown  in  figure  15. 

When  the  aspect  ratio  of  the  [O12]  specimen  is  increased,  the 
thermal  stress  distributions  change.  Figure  21  shows  the  distributions 
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Figure  21  -  Nondimensional  Normal  Stresses  o\^/oj  and  a 
[0,,]  Laminate;  Specimen  Aspect  Ratio  Equa' 


of  and  cTy  stresses  due  to  a  constant  AT  for  the  [O12]  specimen  with 
an  aspect  ratio  of  4. 

Figure  22  shows  the  distribution  of  x=a/2  plotted  against 

y/a  for  a  [O12]  specimen  with  an  aspect  ratio  of  10.  This  figure  shows 
that  the  stresses  have  not  yet  reached  a  uniform  value  at  the  center 
region  of  the  plate.  This  fact  is  further  illustrated  1n  figure  23, 
where  is  shown  plotted  against  x/a  for  a  [O12]  specimen  with 

aspect  ratios  of  2,  4,  6  and  10.  The  stress  is  not  shown  in  these 
latter  two  figures  because  it  is  relatively  small  in  comparison  to  o^. 

Figure  24  displays  the  distribution  of  the  normal  stress  near 
the  free  edge  of  a  [O12]  specimen  with  an  aspect  ratio  of  10.  This 
distribution  is  similar  to  that  seen  previously  for  the  [9O12]  laminate. 
However,  since  for  the  [O12]  laminate  is  a  positive  number  for  a 
heated  specimen,  tensile  values  of  the  a  stresses  will  develop  near 

A 

the  corners  when  the  [O12]  specimen  is  heated.  It  again  should  be  noted 
that  the  value  of  used  to  nondimensional ize  the  previous  results  is 
different  for  the  two  laminates.  The  ratio  of  for  the  [O12]  laminate 
to  the  value  of  0^  for  the  [90.J2]  laminate  is  -9.52.  Taking  this  ratio 
into  account,  the  values  of  the  a  stresses  that  appear  in  the  corners 
of  both  the  specimens  are  of  the  same  order  of  magnitude,  but  have 
different  signs.  When  the  specimens  are  heated,  tensile  stresses  appear 
near  the  corners  in  the  [0.|2]  laminate  while  compressive  stresses  appear 
in  the  [90^2^  laminate. 

When  the  [0-j2]  laminate  is  heated,  the  specimen  must  resist 
tensile  stresses,  stresses  that  act  in  a  direction  perpendicular  to  the 
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fiber  direction;  it  would  appear  that  the  [90^2^  specimen  is  a  more 
satisfactory  test  configuration  than  the  [0-123  laminate. 

Finally,  figures  25  and  26  illustrate  the  typical  laminate 
deflection  behavior  due  to  uniform  heating  for  [0^2]  [9O12] 
specimens  with  aspect  ratios  of  4. 


Figure  25  -  Specimen  Deflection  Figure  26  -  Specimen  Deflection _ 

Pattern  Caused  by  Heating  to  Uniform  Pattern  Caused  by  Heating  to  Uniform 
Temperature.  [0-123  Laminate;  Specimen  Temperature.  [90i23  Laminate;  Specimen 

Aspect  Ratio  Equals  4.  Aspect  Ratio  Equals  4. 
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5.0  CONCLUSIONS 


An  analytical  study  has  been  presented  that  provides  information 
about  the  stresses  and  deflections  caused  by  mechanical  loading  and 
thermal  expansion  of  a  linear  elastic,  laminated  composite,  rail  shear 
test  specimen.  Results  indicate  that  rail  shear  testing  of  the  [9O12] 
graphi te/polyimide  laminate  should  provide  accurate  experimental  data 
for  the  determination  of  the  shear  modulus  of  the  laminate,  even  for 
specimen  length-to-width  ratios  as  small  as  3.  Significant  experimental 
error  is  to  be  expected  for  other  symmetrical  laminates  such  as  the 
[±462]^,  the  [90/+45/0]  laminate  or  the  [O12]  laminate,  unless  length- 
to-width  ratios  larger  than  6  are  used. 

Regions  with  large  shear  stresses  and  normal  stresses  are  found  to 
occur  near  the  corners  of  the  specimens  when  they  are  subjected  either 
to  mechanical  loading  or  to  uniform  heating  or  cooling.  This  indicates 
that  the  rail  shear  test  may  not  be  entirely  satisfactory  for  determining 
laminate  allowable  shear  stress. 

Analysis  of  the  effect  of  heating  an  assembled  rail  shear  specimen 
to  a  uniform  temperature  shows  that  the  [O12]  laminate  is  more  adversely 
affected  by  differential  expansion  effects  than  is  the  [9O12]  laminate. 

In  both  cases,  regions  of  large  normal  stress  appear  near  the  corners. 
However,  for  large  aspect  ratio  specimens,  these  stresses  are  less 
severe  in  the  case  of  the  [9O12]  laminate  than  they  are  for  the  [O12] 
laminate. 
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APPENDIX 


LAMINA  GEOMETRY  AND  MATERIAL  PROPERTIES  USED  IN  THE  PRESENT  STUDY 

A.  Graphite/Polyimide 

E^  =  149.6  GPa  (21 .7  x  lO^psi) 

=  8.273  GPa  (1.2  x  lO^psi) 

G^2  =  4-137  GPa  (0.6  x  lO^psi) 

^12  =  0.27 
^21  ”  0-015 
a-|  =  0 

a2  =  26.1  X  10"^  m/in/°C  (14.5  in./in./op) 
a  =  7.62  mm.  (0.30  inches) 
t  =  0.1016  mm/ply  (0.004  inches/ply) 

B.  Titanium 

E  =  113.8  GPa  (16.5  x  10^  psi) 
y  =  0.342 

a  =  9  X  10"®  m/m/°C  ,(5  x  10"^  in./in./°F) 

Rail  width  =  2.54  cm.  (1  inch) 

Rail  depth  =  2.54  cm.  (1  inch) 
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